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Abstract 

Using a macroscopic analysis, we demonstrate that time-dependent noncollinear spin transport 
may show a wavelike character. This leads to modifications of pure spin-diffusion dynamics and 
allows one to extract a finite spin-signal propagation velocity. We numerically study the dynamics 
of a pure spin current pumped into a nonmagnetic layer for precession frequencies ranging from 
GHz to THz. 
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I. INTRODUCTION 



Transporting information encoded in electronic spins through layers of ferromagnetic and 
normal metals is a central theme of magnetoelectronics.- Structures, in which all spins are 
are essentially collinear, i.e., parallel or antiparallel, have been thoroughly investigated in 
experimental and theoretical studies. The quasi-static properties for the special case of 
structures with collinear spin and magnetization directions where the spin-polarized cur- 
rent flows perpendicularly to the plane of the layers,- can be analyzed in terms of a scalar 
space-dependent spin accumulation for up and down spins.-'- The functionality of collinear 
magnetoresistive structures can be enhanced by including tunneling elements.-i^>^ Although 
collinear spin transport is of importance for certain variants of giant and tunneling magne- 
toresistance effects, a non- collinear alignment of spin and magnetization orientations leads 
to additional degrees of freedom for the manipulation of spin angular momentum and has 
attracted much attention in recent years.-^'^ For instance, one can exploit the angular de- 
pendence of the giant magnetoresistance effect^, or can change the the alignment of spins 
by spin currents, leading to the phenomenon of spin transfer torque^'i^'i^'i^ and potential 
novel applications.—. A different method to exploit the freedom of noncoUinear spin ori- 
entations in magnetic multilayers is the use of magnetization precession in a ferromagnetic 
layer, which "pumps" a spin currents into an adjacent nonmagnetic metal. A precessing 
magnetization, which is necessary for spin pumping, creates the need to deal with a time- 
dependent orientation of the spins in the whole multilayer, so that it becomes essential to 
study dynamical noncollinear spin transport problems. 

We are concerned with a theoretical analysis of the propagation of signals encoded in 
a spin current, which flows through a multilayer structure with noncollinear magnetization 
and spin directions. Most investigations of time-dependent noncollinear spin transport are 
based on the Bloch-Torrey diffusion equations for the nonequilibrium magnetization or spin 
accumulation.^^ These equations essentially describe spin transport as a diffusion process and 
therefore show the same problem as the spin diffusion equation for collinear spina^'^^'i^i^: 
no finite propagation velocity for a spin signal can be defined because the diffusion equation 
leads to a finite spin current density everywhere as soon as there is a source. Recently, we 
showed that this difficulty can be resolved for collinear spin transport by using a "telegraph" 
equation, which generalizes the diffusion equation, and leads to noticeable differences from 
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the diffusion equation results for frequencies exceeding several 100 GHz for metals such as 
copper.— Importantly, the telegraph equation shows a wave-diffusion duality, which enables 
one to define a finite propagation velocity for the spin signal. In this paper, we use a similar 
treatment for noncollinear spin transport to show how a finite signal propagation velocity 
arises in this case. We predict that noncollinear spin transport at high frequencies shows a 
dynamics that is more complicated than what is expected from an analysis using the spin- 
diffusion equation. We numerically analyze the propagation of a spin current pumped into 
a nonmagnetic metal by a precessing magnetization in an adjacent ferromagnetic layer. 

This paper is organized as follows. In Sec. [Tll we present the macroscopic dynamical 
equations governing noncollinear spin transport. In Sec. Illlt the dynamical equations are 
combined into a telegraph equation, which is studied analytically to discuss qualitative 
aspects of dynamical noncollinear spin-transport. In Sec. llVt we solve numerically the 
dynamical equations for the spin transport, and the main conclusions are summarized in 
SeclVl 



II. DYNAMICAL EQUATIONS 

In nonmagnetic conductors and some ferromagnetic metals,— the dynamics of conduction 
electrons under the influence of external fields can be described by a generalized semiclassical 
Boltzmann equation^^"^ 



,df) i { de dp\ i \ de df)\ ,dp 



dt dr\ 2\dr' qI\ ^ ' dt 



(1) 



col 



which we take as the starting point for our analysis of time-dependent noncollinear electron- 
spin transport in these systems. In Eq. ([1]), p(r, k^t) is the single particle density matrix in 
spin space, 

P=f'""'"V (2) 

\Pi\ Pa j 

e{f, k,t) is the effective single-particle energy matrix, and {■, ■} and [■, ■] denote respectively 
the anticommutator and commutator for matrices in spin space. For completeness, we 
remark that in Eq. ([2]), the single-particle density matrix 
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is defined by a statistical average over creation and annihilation operators and c, with 
normalization volume V. The diagonal matrix elements and pn are the electron distribu- 
tion functions of the spin-up and spin-down, respectively, whereas the off-diagonal elements 
Pu ~ Pit represent the spin coherence.-- Because the unit matrix / and the Pauli matrices 
a^c, d'y, dz form a basis for 2x2 matrices, the spin-density matrix p can be represented by 
P= (l/2)[(PTT+Pii)^ + ^-o-], where M = Tr(6-p) = (2Rep||, -2Imp||, -p||) is the Bloch 
vector and <t the vector of Pauli matrices. 

Before proceeding from Eq. ([1]) for the spin-density matrix to equations for macroscopic 
quantities, such as spin currents and spin accumulation, we list a few assumptions made 
about quantities occurring in Eq. ([T]). First, we consider only layered structures whose 
extensions perpendicular to the growth direction [x axis) are infinite, and we also assume 
that the electric fields E = Ex/\x\ is oriented along the growth direction x. Second, the 
effect of magnetic fields on the orbital motion of electrons is neglected. These magnetic fields 
include the static external magnetic field Bs and the magnetic field generated by induction 
due to the time- dependent electric field E{x,t).— We therefore assume that the electric 
field E{x, t) = —d(f){x, t)/dx can be derived from a time-dependent electric potential t). 
Third, an isotropic effective mass model for the spin-degenerate condutcion electrons is used, 
i.e., Ek = h'^k'^/{2m*) = m*v'^/2, where k and v denote the the electron wave vector and 
velocity, respectively. Thus we have to deal with a spin density matrix p that depends only 
on X and has cylindrical symmetry around the x axis in k space. 

Finally, we make a relaxation-time approximation for the collision term^ 

dp 

col ^ ^1 



P - {P)a {P)a - {i/2)TT{pU^ 



dt 

where r and Ti are the momentum and spin relaxation times, respectively. Moreover, 
(p)a = (47r)^^ J dQj:p is the angular average in the momentum space. By using Eq. (jlj) for 
the collision term, we have assumed that the longitudinal spin relaxation time Ti is equal to 
the transverse one T2. The validity of this approximation is discussed in detail by Ref. 
Note that Ti in Eq. (jlj) is one half of Tgf used in Eq. (2) of Ref. [27. 



With above simplifications, the effective single-particle energy e{f, k, t) is simplified to 

i{x, \v\,t) = sqI + is, where Eq = h'^k'^ /{2m*) — e(j){x,t) and is = —i-i ■ Bg = Pbct ■ Bs. 

Therefore, Eq. ([T]) simplifies to 

dp dp eEdp I p - (p), (p). - (//2)Tr(p), 

^ + ^— + -7(m X 5,) ■ cr = , (5) 

at ox m* ovx 2 r Ti 



where 7 = gfis/h is the absolute value of the electron {g ^ 2) gyromagnetic ratio. 

To derive macroscopic spin transport equations comparable with the Bloch-Torrey diffu- 
sion equation, we need to sum over the electron wave vector k or, equivalently, the velocity 
V in Eq. We first derive an equation for the spin density^^'^^ by multiplying both sides 
of Eq. by o'/V, taking the trace, and summing over v 

dns{x,t) - / ,N 5 ns{x,t) djs{x,t) 

^^ = -7n.(x,t)x5.-^^ (6) 

where n^(x,t) = V^"^ Et;Tr(crp) = X],;^* and /.(x, t) = 1/"^ X]^^xTr(o-p) = V-'^Y.v'^xU 
are the spin density and spin current density, respectively. For the spin current density, 
we multiply both sides of Eq. ([5]) by v^o-fV, take the trace, and sum over v. Using the 
expansion ( ]A2|) for the velocity dependence of the spin density matrix and the procedure in 
Appendix A, we obtain 

Ux, t) = -D ^'^'gl'^^ - f^E{x, t)ns{x, t) - T^U^, t)xBs- r^^£^, (7) 



where 



D = fr (8) 



is the diffusion constant and fi = eT/m* the electron mobility. Note that ns{x,t) and 
js{x,t) defined above are the particle (electron) number densities, which can be converted 
to the charge, spin, and magnetic moment densities by multiplication with — e, h/2, and 
—/is, respectively. The spin density ns{x,t) can also be converted to the chemical potential 
difference fis{x,t), i.e., the spin accumulation, by the relation ns{x,t) = M^isix^t), where 
M = 47rm*^f / is the density of states at the Fermi level of the electron gas for one spin 
orientation.— 

Equation ([7]) resembles the dynamical equation for the spin current derived by Qi and 
Zhang2^ using a "mean field" approximation. Our derivation shows that their quantity 
is equal to Vp/3. As will be discussed in the next section, this is the wavefront velocity for 
a spin disturbance, which plays an important role in spin-signal propagation dynamics—. 

III. TELEGRAPH EQUATION 

To see the physical significance of Eqs. ([6]) and ([7]) for the time- dependent noncollinear 
spin transport and compare them with the Bloch-Torrey equation, we combine them by 
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eliminating js{x,t) into a form reminiscent of a telegraph equation^i 



dt^ \ Ti' dt tT, ^ 



d A 1 , 

2a + <7 + iT' 



ns{x,t) X Bs 



+ 72[n,(x,t) X Bs] X Bs 
+ 



(9) 



2d'^ns{x,t) ^ fiE{x,t) dns{x,t) ^ fidE{x,t) 



+ 



ns{x,t). 



* dx^ T dx T dx 

Similarly, one can also derive a telegraph equation for js{x,t) by eliminating ns{x,t) from 
Eqs. ([6]) and ([7]). Equation (Q contains a second-order time derivative, which is absent in 
the spin diffusion equation. The second-order time and space derivatives lead to a wave 
character in addition to its diffusion character, and thus yield a well-defined propagation 
velocity Cs for the signal in time-dependent noncollinear spin transport in a similar way to 
the collinear case.— 

Assuming the static magnetic field Bs to be oriented along the z axis and separating the 
components perpendicular (transverse) and parallel (longitudinal) to Bg in Eq. (Q, we have 



d n 



+ 



1 1 

T Ti 



dn, 



dt 



+ 



n 



x{y) 



+ i-hBs 



\9 [I 1 



dx"^ 



+ 



fiE dn 



x{y) 



T 



+ i^^n^iy) 
dx T dx ^ 



and 



d'^nl 



1 1 

T Ti 



dnl 



.d'^nl 



+ 



^E dnl 



II dE , 

T) 

Tdx 



(10) 



(11) 



dt tTi ' dx^ ' T dx 
In the following, only the equation for the transverse component [Eq. 01 01) ] will be discussed, 

since the equation for the longitudinal component is similar to that of the collinear case.— 

For vanishing electric field, i.e., E' = 0, we seek damped and dispersive wave solutions to 

Eq. of the form 



n^{x,t) = nQexp[i{kx — ut)], 
n^(x,t) = n()exp[i{kx — ut + 



(12) 
(13) 



where u is the angular frequency and k = kr + iki the complex wave vector. Substituting 
Eqs. f lT2|) and f|T3|) into Eqs. ffTOj) . we obtain the dispersion relation 

+ tuj{l/T + 1/Ti) - l/(rTi) - clk^ + -f^Bl - 7^, [2u + i (1/r + 1/Ti)] sin0 = 0, (14) 

where is restricted to = ±(7r/2) + 2ri7r and n is an integer, because and n\ must 
satisfy the system of equations (fTOjl at the same time. According to Eqs. ([T2]) and (fT3j) . 



6 



= +(— )7r/2 corresponds to the rotation direction of the transverse component of ns{x,t) 
with X at time t. For definiteness, we study the case with (p = 7i/2 in the following. 
Substituting k = kr + iki into Eq. (fn|) and separating the real and imaginary parts, we have 

1 



/,2 



+ ul^a^ + {-)b 



(15) 



where = u — -^Bg and h = cj^g — ^. Here, the constants a = 1/t + 1/Ti and ^ = l/(rTi) 
have been introduced. The wavelength and damping length can be defined as A = 211 /k^ 
and Id = l/ki, respectively. The equation of the critical angular frequency cUcrit, above which 
the wave character is significant, can be derived by setting A = Id, 

^ Uil + ri) + v/52(l + r/)2 + 4r/l ^ 5 + {6 + h^, (16) 



, ,crit _ 
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where S = n — l/(47r) ^ 3.06 and r] = r/Ti. Then, we have UcntT = 3.06 + 3.4// + t'jBs 
approximately. 

IV. DYNAMICS OF PUMPED SPIN CURRENT 

In this section, we study the evolution of the spin current injected into a nonmagnetic 
layer by the spin-pumping mechanism.-l^ In a junction composed of a ferromagnetic {x < 0) 
and a nonmagnetic {x > 0) layer, the magnetization precession of the ferromagnet around an 
external magnetic field -Bpump acts as a "spin pump" which transfers spin angular momentum 
from the ferromagnet to the adjacent nonmagnetic layer. The spin current density pumped 
into the nonmagnetic layer igi^i^EiSS 

- tS- - f • ('^' 

where g^^ is the spin-mixing conductance and 5* the area of the interface. Here, rh is the 
unit vector for the magnetization of the ferromagnet. Note that the pumped spin current 
has been converted to a particle number current density j^"™^. Since we are interested in the 
spin current pumped into a nonmagnetic layer and not in the dynamics of the ferromagnet, 
we neglect the back-flow spin current f^^^^^ which flows from the nonmagnetic layer to the 
ferromagnet due to the spin accumulation in the nonmagnetic layer.29 Although the back- 
flow spin current can limit the achievable spin current into the nonmagnetic conductor, we 
do not approach this limit here. With this simplification, we have JP"™? = J^(a; = 0,t), 



where js{x = 0,t) is the spin current density at the left boundary of the nonmagnetic layer. 
Separating the components perpendicular and parallel to the magnetic field -Bpump, we can 
write js{x = 0, t) as 

j^{x = 0,t) = g^^{4nS)-^uJsm{2e)cos{ujt) (18) 
jy{x = 0,t) = g^'^{AnS)-^usm{2e)sm{ut) (19) 
f^{x = 0,t) = g^^{2'KS)-^uJsm^e, (20) 

where u is the angular frequency of both the magnetization precession and the spin current 
density js{x = 0,t). Here, ut is the angle between jf (j^ and j^) and the x-axis. 9 is the 
angle between m and -Bpump, and meanwhile 9 is also the angle between js{x = 0,t) and 
xy-plane. The amplitude of is much larger than j^, since 9 is very small under the usual 
radio-frequency excitation conditions.— Therefore, we will focus on jj- in the following. 

The propagation of jf{x = 0,t) into the nonmagnetic layer is described by Eqs. (jH]) 
and (171). In a typical setup for spin pumping, there is no electric or magnetic field in the 
nonmagnetic layer, i.e., E = and Bg = 0. Now, separating the components perpendicular 
and parallel to the magnetic field -Bpump, we can rewrite Eqs. IQ and ([7]) as 

dnt dj+ _ _nt . . 

at ^ dx ~ Ti ' ^^^^ 

where = + and j+ = + ijf are introduced to simplify the notations. The 
equations for the parallel component can be obtained after replacing and j'^ by and 
jg in Eqs. ( pTl) and ( 122|) . respectively. The method of characteristics used for the numerical 
solution to Eqs. (I2T1) and (l22l) is outlined in Appendix [B] 

In our numerical calculation, Cu and permalloy (Py) are chosen as the materials for 
the nonmagnetic and ferromagnetic layers, respectively. The Fermi velocity of Cu is vp = 
1570nm/ps and thus the wave-front velocity is Cs = vp/y/S = 906nm/ps. The momentum 
and spin relaxation times are r = 0.07ps and Ti = 3.5 ps, respectively. The critical frequency 
can be estimated to be i^crit = t^crit/(27r) = 7.11 THz from Eq. (|T6l) . We study several 
pumping frequencies: Ua = I/Tq = 2 GHz, Uf, = l/Tb = 20 GHz, Uc = 1/Tc = 200 GHz, 
and = '^/Td = 8.33 THz. For a Py/Cu junction, g^^S~^ is on the order of 10^^ cm~^. 
The precession cone angle 9 can reach 15° for a sufficiently intense radio-frequency field.— 
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FIG. 1: Snapshots of the spin current density at t = 5/4 Ta, 5/4 T^, and 5/4 Tc, for the 
frequencies, Ua, ft, and fc, respectively (see text), is plotted as vector starting from its x- 
coordinate. 

Therefore, we choose the amplitude of J^, i.e., g^^{4TTS)^^uj sm(26), to be 5 x 10^^nm~^ps~^ 
for the frequencies mentioned above. 

Figure [1] shows snapshots of the spin current density jj- at t = (5/4) T^, (5/4) Tf,, and 
(5/4) Tc, for the frequencies, Ua, ^b, and Uc, respectively. According to Eqs. f|T8|) . jj-{x = 0,t) 
points in the direction of the y-axis at t = 5/4:Ta(b,c), which can also be seen in Fig. ([1]). 
Figure [T] (a) shows that points along y-axis nearly at all x points except that it deviates 
from the y-axis slightly at positions far away from x = 0. The results in Fig. [H (a) are 



approximately consistent with those obtained from the diffusion equation in Refs. 



29 



and 



30 



where it is shown that both the spin current and spin accumulation point along the same 
direction at all positions for all frequencies at certain time point t. This agreement means 
that the diffusion equation provides a good description in the low frequency range.— The 
deviation of from the y-axis at x > increases with frequency and becomes noticeable 
at z/fc = l/Tfc = 20 GHz as shown in Fig. [T] (b). Therefore, the applicability of the diffusion 
equation is questionable in this frequency region. At even higher frequency, Uc = 1/Tc = 
200 GHz, the deviation becomes significant and the diffusion equation is not applicable. 
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FIG. 2: Snapshots of the spin density for the same parameters as in Fig. [TJ 

Moreover, the damping length of jj- decreases with frequency due to the 'skin' effect.— We 
can conclude that the spin diffusion equation is applicable only in the low frequency range 
and amounts to an adiabatic approximation: the external perturbation is assumed to be 
much slower than the internal dynamics of the electronic system. 

Figure [2] shows snapshots of the spin density nj- for the same parameters as in Fig. [H 
The spin density nj- deviates from y-axis at x = and is noncollinear with at x > at all 
of the three frequencies. This feature is different from the result of the diffusion equation, 
where jj- and are collinear.— The phase shift and the amplitude of also vary with 
frequency. Moreover, the damping length of decreases with frequency again due to the 
'skin' effect. 

According to Eq. (fTGl) . the diffusion character is dominant at the frequencies considered 
so far, because they are still much smaller than the critical frequency z^crit- This conclusion 
is supported by the numerical results presented in Figs. [T] and [21 although Figs. [T] (c) and [2] 
(c) have already shown weak wavelike character. The deviation from the diffusion equation 
depends largely on the frequency of the spin signal and momentum relaxation time, which 
varies with material, temperature, doping and excitation condition. In the following, we 
show the numerical results for a frequency Ud = 8.33 THz, where the wave character is 
significant according to Eq. (ITB]) . 
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FIG. 3: Snapshots of at t = T^, 2Tii., and ST^, where Td = 0.12 ps. The sohd (dashed) curve is 
for j| (jf). 

Figure [3] shows snapshots of the spin current density aX t = ITd, 2Td, and ST^, re- 
spectively. The wave form and wave front are clearly visible in Fig. [31 The propagation 
velocity of the spin signal can be estimated by tracking the motion of the wave front. The 
result is approximately equal to the analytical result Cg = 906nm/ps. The phase velocity 
can also be estimated by measuring the wavelength A and using Vp = X/Td. The result 
is roughly equal to the wave front velocity Cg, which also indicates the significance of the 
wave character, albeit on the length scale of the damping length (dynamical spin diffusion 
length). To demonstrate the wave character more directly, we plot the results of Fig. [3] (a) 
again in Fig. |U where jj- is shown in a vector plot. Note that Ud is beyond the frequency 
range in which Eq. (fTTI) is valid, because Eq. (ITTI) is only applicable in the adiabatic limit, 
u ^ 1/r.— Unfortunately, there is no corresponding theoretical result for the nonadiabatic 
spin-pumping in the literature. However, it is a reasonable guess that the pumped spin 
current density in the nonadiabatic regime preserves the basic feature of Eq. (fTTj) : jf"™^ 
rotates with a certain fixed frequency. Therefore, the spin current predicted by our results 
should be at least qualitatively accurate in this frequency range. 
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FIG. 4: Snapshot of (plotted as vector) at t = 1 T^. 

V. SUMMARY 

We showed that time- dependent noncolhnear spin transport exhibits a wave character for 
modulation of the spin current on timescales shorter than an inverse critical frequency. A 
finite propagation velocity for the spin signal can be defined due to this wave character. The 
spin diffusion equation is recovered only for modulation with frequencies less than the critical 
fequency, and amounts to an adiabatic approximation of time-dependent spin transport. 
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APPENDIX A: DERIVATION 



Equations (6) and (7) of Ref. l27| are derived using the "mean field" approximation 



'^vl{dp/dx)^vlY,{dp/dx). (Al) 

V V 

Here we show that = by evaluating the sums occurring in Eq. (jXTD. We start with 
the LHS, which we denote by Ji = Ylv'^xi^p/dx). Due to the cylindrical symmetry of the 
system around the x axis in velocity space, p can be expanded in Legendre polynomials of 
u = cos 6*, where 6 is the angle between v and the OC dbXlS^ clS 

oo 
n=0 
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Transforming the summation into an integral, we have 

h = — j-^ — / duu / dvv 2_^—pn[v,x)Pn[u). (A3) 
J ~i Jo 

Using = [2P2{u) + Po(w)]/3, we write the integral as 

/i = ^^/ / t/«-[2P2H+Po(«)]PnH. (A4) 



Making use of the orthogonality relation of Legendre polynomials, we have 



^ . L dx 



4 2 

— P2{V,X) + -po{v,x) 



(A5) 



If the system is weakly anisotropic, we can neglect the second-order term P2{v, x) 



h^^^ dvv-po{v,x). (A6) 



This approximation is consistent with Ref. 



4|, where the second-order term of the Legendre 



polynomials is neglected and it is shown that this is valid if -^/r/ (2Ti) ^ 1. 

Because dpo{v,x)/dx is zero unless v falls in a small region [vp — Av,vf + Av] around 
the Fermi velocity fi? of a system with a degenerate electron gas, we have approximately 

/i = — / dvv -w-Po{v,x) = — — —— / dvv Tr:Po{v,x) . (A7) 



3 JvF~Av dx ' 3 Jq dx 

We now need to evaluate the RHS of Eq. (lAip , which we denote by 

l2 = vl—^2nj du J dvv^—p{v,x)=v^ — — — J dvv^—po{v,x). (A8) 

where, in the last line, we used that the integral over u projects the contribution of Pq out 
of p{v,x). Because Ji = I2, we conclude that = Vp/3 = c^. 



APPENDIX B: NUMERICAL SOLUTION 



The basics of our numerical method have been outlined in Appendix A4 of Ref. l21 
For present calculation, it has to be augmented by a discretized version of the boundary 
condition on at ferromagnet /nonmagnet interface, 

{At/T, + 2) = - {At/T, - 2) n+i, + c;^ {At/r - 2) jtj^, + c;^ ( At/r + 2) jtt' 

(Bl) 
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where the subscripts i and superscripts / stand for the discrete space-time points, and At is 
the numerical time step. 



Electronic address: yaohuizhu@gmail.coin 



^ URL: http://www.physik.uiii-kl.de/sclmeider 



1 A. Brataas, G. E. W. Bauer, and P. J. Kelly, Phys. Rep. 427, 157 (2006). 

^ W. P. Pratt, Jr., S.-F. Lee, J. M. Slaughter, R. Loloee, P. A. Schroeder, and J. Bass, Phys. Rev. 

Lett. 66, 3060 (1991). 
3 P. C. van Son, H. van Kempen, and P. Wyder, Phys. Rev. Lett. 58, 2271 (1987). 
^ T. Valet and A. Pert, Phys. Rev. B 48, 7099 (1993). 

^ G. Schmidt, D. Ferrand, L. W. Molenkamp, A. T. Fhip, and B. J. van Wees, Phys. Rev. B 62, 
R4790 (2000). 

6 E. I. Rashba, Phys. Rev. B 62, R16267 (2000). 

7 A. Pert and H. Jaffres, Phys. Rev. B 64, 184420 (2001). 

^ Y. Tserkovnyak, A. Brataas, G. E. W. Bauer, and B. I. Halperin, Rev. Mod. Phys. 77, 1375 
(2005). 

9 J. Barnas, O. Baksalary, and A. Pert, Phys. Rev. B 56, 6079 (1997). 
^° J. C. Slonczewski, J. Magn. Magn. Mater. 159, LI (1996). 
" L. Berger, Phys. Rev. B 54, 9353 (1996). 

12 M. Tsoi, A. G. M. Jansen, J. Bass, W.-C. Chiang, M. Seek, V. Tsoi, and P. Wyder, Phys. Rev. 
Lett. 80, 4281 (1998). 

13 E. B. Myers, D. C. Ralph, J. A. Katine, R. N. Louie, and R. A. Buhrman, Science 285, 867 
(1999). 

1^ J. A. Katine and E. E. Fullerton, J. Magn. Magn. Mater. 320, 1217 (2007). 
1^ Y. Tserkovnyak, A. Brataas, and G. E. W. Bauer, Phys. Rev. Lett. 88, 117601 (2002). 
1^ L Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 76, 323 (2004), in particular. Sec. 
IILA.l. 

S. Zhang and P. M. Levy, Phys. Rev. B 65, 052409 (2002). 

18 E. I. Rashba, Appl. Phys. Lett. 80, 2329 (2002). 

19 J. Zhang and P. M. Levy, Phys. Rev. B 71, 184417 (2005). 

14 



L. Cywinski, H. Dery, and L. J. Sham, Appl. Phys. Lett. 89, 042105 (2006). 

Y.-H. Zhu, B. Hillebrands, and H. C. Schneider, Phys. Rev. B 78, 054429 (2008). 

J. Zhang, P. M. Levy, S. Zhang, and V. Antropov, Phys. Rev. Lett. 93, 256602 (2004). 

H. Smith and H. H. Jensen, Transport Phenomena (Clarendon Press, Oxford, 1989). 

J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

L. J. Sham, J. Magn. Magn. Mater. 200, 219 (1999). 

N. W. Ashcroft and N. D. Mermin, Solid State Physics (Brooks/Cole, Belmont, CA, 1976). 

Y. Qi and S. Zhang, Phys. Rev. B 67, 052407 (2003). 

A. D. Margulis and V. A. Margulis, Pliysica B 193, 179 (1994). 

Y. Tserkovnyak, A. Brataas, and G. E. W. Bauer, Phys. Rev. B 66, 224403 (2002). 

A. Brataas, Y. Tserkovnyak, G. E. W. Bauer, and B. 1. Halperin, Phys. Rev. B 66, 060404(R) 

(2002). 

T. Gerrits, M. L. Schneider, and T. J. Silva, J. Appl. Phys. 99, 023901 (2006). 



15 



